New exact fronts for the nonlinear diffusion equation with quintic nonlinearities 
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We consider travelling wave solutions of the reaction diffusion equation with quintic nonlinearities 
Ut = u X x + fiu(l — w)(l + ctu + flu 2 + 7U 3 ) . If the parameters a, (3 and 7 obey a special relation, then 
the criterion for the existence of a strong heteroclinic connection can be expressed in terms of two 
of these parameters. If an additional restriction is imposed, explicit front solutions can be obtained. 
The approach used can be extended to polynomials whose highest degree is odd. 



I. INTRODUCTION 

The nonlinear diffusion equation ut = u xx + f(u) models phenomena in diverse fields such as population growth, 
kinetics of phase transitions, chemical reactions and many others. Of special interest is the case when the function 
/ is such that there exist two steady states, one stable and one unstable. We shall assume that the equation has 
been scaled so that the unstable state is u u = and the stable state is u s = 1, and we consider functions / which 
are positive in (0, 1). Then sufficiently localized initial conditions evolve into a travelling front which joins the two 
steady states p]]. The speed at which the front propagates, c* is equal or greater than the linear marginal stability 
value cl = 2\/ /'(0). In many cases the asymptotic speed of propagation is exactly the linear value cl = 1\J '/'(0) 
obtained by the linear marginal stability criteria ^J,[| . There are cases however when the front propagates at a speed 
greater than this value, case which is referred to as that in which a nonlinear speed selection mechanism operates. 
Explicit expressions for this special nonlinear front or strong heteroclinic connection and its speed have been obtained 
for particular choices of /. All the known solutions correspond to functions / of the form f(u) = llu + u n — u 2n_1 
which, for fx positive but smaller than a critical value fj, c , are strongly heteroclinic |?J. The purpose of this article 
is to show, using as an example a quintic polynomial /, that the criterion for the existence of special fronts can 
be formulated in many cases in a simpler way that enables one to decide whether for a certain / there is a strong 
heteroclinic connection even if the exact solution for the front is not known. We find new exact fronts for this quintic 
polynomial for / together with a criterion for strong heteroclinicity in terms of the parameters of the polynomial valid 
even when no explicit solution for the front can be obtained. Similar results can be obtained for polynomials whose 
highest degree is odd. The knowledge of exact solutions is of interest not only as a curiosity, they are also needed in 
the framework of the recent proposal of structural stability || , the knowledge of the speed for a specific form of / 
enables the calculation of the speed for small perturbations to / using renormalization group techniques. 

In section 2 we state the problem and reformulate already known results, and in section 3 we give the new results 
for the quintic polynomial. 

II. MONOTONIC FRONTS OF THE REACTION DIFFUSION EQUATION 

We consider the reaction diffusion equation 

ut = u xx + f(u) 

with f(Q) — 0, /(l) = 0, f'(0) > and / > in (0, 1). Given these conditions on / then there exist fronts that 
connect the unstable fixed point it = to the stable fixed point u — 1. Travelling wave fronts u(x — ct) satisfy the 
ordinary differential equation 

u zz + cu z + f{u) = lim u = 1, lim u = 0, (1) 

z — > — 00 z — >oo 

where z — x — ct and we assume that c is positive. A front joining the stable fixed point 1 to the unstable point is 
monotonic if in addition its derivative du/dz does not change sign. If we search for monotonic fronts it is convenient 
to consider the dependence of z as a function of u, or rather the dependence of v(u) = —(dz/du) -1 as a function of 
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u. For a monotonic solution of equation (1), u(z) decreases monotonically as z goes from — oo to oo, therefore the 
function v(u) is well defined and is positive between 1 and and vanishes at the fixed points. One readily finds that 
the equation for v(u) is 

dv 

v(u)—-cv(u) + f(u) = 0, (2a) 
au 

with 

v(0) = v(l) = 0, with v > 0. (2b) 

Since the endpoints are singular we must determine the behavior near them analytically. If we consider functions / 
analytic around and with /'(0) > 0, then near u = we find 

v(u) = a\u + a 3 / 2 u 3 ^ 2 + a 2 u 2 + a 5 / 2 u 5 ^ 2 + a 3 u 3 + . . . 
where the first terms are given by 

a\ - ecu + /'(0) = (3a) 
5 

«3/2(2 a l - c ) = ( 3b ) 

a 2 (3ai - c) + lf"(0) = (3c) 
7 7 

°5/2(c - -at) - -a 3/2 a 2 = (3d) 

and so on. That the leading term in the expansion of v near zero is linear in u is due to the fact that the front in the 
original coordinates u{z) approaches the fixed points exponentially. Since v must be positive between and f , a\ must 
be real and positive. The two roots for a\ are given by aip = (c+ \J c 2 — 4/'(0) )/2 and a\u = (c— \J c 2 — 4/'(0) )/2. 
The minimum speed c for which there may be a monotonic front is the linear marginal speed value cl = 2-\/ f'(0) 
value at which the roots coincide a\p — a^/ = cl\l- For speeds greater than this value a\M < &il < a ip- Strong 
heteroclinic solutions or special nonlinear front profiles are those associated with a\p. From the expansion at the origin 



it follows from ( pb[ ) that either c = 5ai/2 or a 3 / 2 = 0. In the first case we find that c = 5v//'(0)/6 « 2.041 ^/flfi) 
and ai = v/2/'(0)/3 = aiM- As it is known, these solutions are not a preferred asymptotic state. Strong heteroclinic 
connections can be achieved only if a 3 / 2 = 0, all half integer coefficients vanish then and v(u) = a±u + a 2 u 2 + . . .. 
Near u = 1, assuming /'(l) < 0, 

v(l - u) = 6i(f - u) + b 2 (l - u) 2 + 6 3 (1 - uf + . . . 

where b\ is the positive solution of 

b\ + ch + /'(l) = 0. 

There is only one positive solution for b\ , the rest of the coefficients follow easily. 

It is convenient to introduce a new parameter A defined by c = \a\. It is not difficult to realize that whenever 
1 < A < 2 then the solution for v is strongly heteroclinic, that is, associated with a\p and when A > 2 it becomes 
associated with q\m] hence for A > 2 the linear marginal speed is selected. If c = \a\ then 



\ //'(°) //'(°) U s 

^ a - !' fll = VA-T- (4) 



At A = 2, the speed attains its linear value cl — 1\J '/'(0). The problem then is to determine the value of A. This 
transition value A — 2 is not associated with any specific nonlinearity, it is valid for any / which satisfies the conditions 
given above. 

All exact front solutions given in the literature, P,p|,pI)Pj7|] correspond to functions / for which an exact solution 
for v is of the form 
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v n (u) = aiu(l — u n 1 ) 

which is an exact solution of equation @ for 

/«(«) = /'(o) (« + (1 +^ A) ^ - x^" 2 - 1 ) (5) 

We observe that for A = n + 1 we recover the solutions of Kaliappan since n > 1, A is greater than 2, so none of 
them are strongly heteroclinic. The front corresponding to v n is given implicitly by 

''" (6) 



v n (u 



and explicitly by @,|TT| 

u n (z) = 



(1 -f c -(n-l)zai);r=T 



The criterion for the existence of strongly heteroclinic fronts together with their exact expression has been given 
[Q] for functions / of the form f(u) = \xu + u n — u 2n ~ l . The critical value for \x given in (?J for the transition from 
a strong heteroclinic connection to a simple nongeneric connection (a solution associated with cjim) is equivalent to 
the value A = 2 after suitable rescaling. It is perhaps convenient to see it in the example given by Van Saarloos 



k = 0xx + <j> + d(/) 3 - 4> b 



which has a strongly heteroclinic connection for d > 2/v3, of speed 



t -2 + 2 % /4Td I 
v = f • 

To identify the value of A from equation (||) we must scale the equation for <fi so that the stable point is at 1. To 
do so we let u = K<fr, u satisfies u t = u xx + u + K 2 du 3 — K 4 u 5 and the stable state is u = 1 if 



K 



, d + Vl+d 2 



where the positive sign is chosen to obtain a real K. Now we compare with equation (||) with n = 3. We see that 
/'(0) = 1 and that 



. , 3d 

A = 4 - 

K 2 



or, in terms of d, 



_ 4V4 + d 2 -2d 

A — 



d+V^+d 2 

It is straightforward to see that the critical value d = 2/V3 is exactly A = 2 and that the speed 

A 

c = , 

x/A^T 



III. NEW SOLUTIONS 

Now consider fronts for / being a quintic polynomial in u. This problem was considered in JTcfl but no explicit 
solutions were found and no attempt to examine the conditions for the transition from the linear to the nonlinear 
regime were made. Here we show under which conditions a closed form can be obtained, together with some examples 
and the condition for strong heteroclinicity A < 2 in terms of the parameters of the function /. Evidently the value 
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of the parameter A can be determined analytically only if an exact solution for v is known. The most general form of 
a quintic polynomial that vanishes at and 1 is 

f(x) = fix(l - x)(l + ax + (3x 2 + -fx 3 ) (7) 

where fx, a, (3 and 7 are four arbitrary parameters whose only restriction is given by the requirement /'(0) > and 
/ > in (0, 1). On the other hand, the most general closed form solution for v given a quintic / is given by 

v(u) = a\u(l - u)(l + bu) (8) 

where 6 > — 1. Introducing again the parameter A given above, so that a\ and c are given by equation (|]), equation 
(0) is the exact solution of equation (g) with 

*/ \ titf* m \ fi , (2 + A6-36) 6(5-26) 2 36 2 4 \ 

/(U) = /'(0)u(l - U) ( 1 + ^ — 11 + ^_ 1 V + y-jU 4 J . (9) 

In the solution for v we have three adjustable parameters, A, 6, f'(0) whereas in the most general form for /, four 
adjustable parameters exist. Hence, an exact solution for v can be found chosing three parameters of / arbitrarily 
and the fourth one in terms of them. Choosing /it, and 7 arbitrarily, we identify 

/'(0) = M 



(3/3 + 2 7 )2 

and 

57 



3/3 + 2 7 

and the exact solution exists if 

(2 + Xb - 36) 

a = 

A - 1 

For any other value of a a closed form solution does not exist and we cannot determine the value of A. The criterion 
for the solution to be strongly heteroclinic 1 < A < 2 is expressed now in terms of the free parameters (3 and 7. 

Now we show that an explicit solution for the front in the original coordinates exists only if an additional condition 
on 6, hence a relation between the free parameters /3 and 7 is satisfied. Proceeding as above in equation (^) we find 
that u(z) is the solution of 

,,1+6 

-(6+1) O! z = » (1Q) 

Writing b — n/p the equation for u is 

n n +P 

e -(n+p)a lZ _ ^ 

(1 -u)P{l + bu) n ' K ' 

This can be inverted to obtain the explicit solution for u(z) if n +p = 2, 3, 4. The detailed inversion of all the solvable 
cases is not instructive, here we give one example. Choose n = 2,p = 1, then 6 = 2, and the front is a solution of the 
cubic equation 

u 3 (l + 4e~ 3ai2 ) - 3ue~ 3aiZ - e~ 3cilZ = 0. (12) 
This cubic has two complex roots and a single real positive root which is the desired front, given by 

ni \e^¥- + V4 + e 3aiZ J 3 
u(z) = + ^ j - '— (13) 
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Again this is an exact front for / of the form given by equation (0). It corresponds to a strongly heteroclinic 
connection for A < 2. If one chooses the case n + p = 4 the quartic equation that arises has a pair of complex 
conjugate solutions, a negative solution and a positive solution which is the desired front. For values of b which do 
not allow the obtention of the explicit form of the front u(z) we still have the speed selection criteria in terms of the 
free parameters of the polynomial. 

Closed form solutions v(u) for polynomial f's can be obtained only if / is an odd polynomial. In general, if / is a 
polynomial of degree 2k + 1 that vanishes at and 1, there are 2k free parameters (restricted only by the requirement 
of positivity of /), whereas the corresponding closed form solution for v has k + 1 parameters, which implies that a 
closed form for v, and an explicit expression for A is possible if fc — 1 parameters of / are chosen adequately in terms 
of the k + 1 remaining free parameters. 



IV. CONCLUSION 



We have studied the existence of exact strongly heteroclinic fronts for the reaction diffusion with quintic nonlin- 
earities. We find that the use of phase space enables one to characterize the transition from strongly heteroclinic 
to simple nongeneric fronts in terms of a single parameter A which is the ratio between the speed and the rate of 
decay at infinity. The introduction of this parameter gives a unified way in which to describe the type of solutions 
which is independent of the nature of the nonlinearities. The exact value of this parameter cannot be determined 
analytically when the highest nonlinearity is even, if the highest derivatives are odd it can be determined for special 
choices of parameters. In the case studied here, quintic nonlinearities, the value of A can be determined exactly if 
a special relation betwen the parameters of the equation is satisfied. It is not necessary to know the exact solution 
u(x — ct) in order to determine whether a strong heteroclinic connection exists. If an additional restriction on the 
parameters is imposed, new exact solutions can be found. We have illustrated this situation for one particular choice, 
a whole family of exact solutions can be constructed. The use of phase space is not only useful as an aid to find exact 
solutions, it can be used to obtain a lower bound on the speed, valid for all /, which allows one to determine the 
range of parameters for which strongly heteroclinic connections exist |12j . 
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